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We have found ultraviolet asymptotic slutions of the Schwinger-Dyson equation for the gluon and ghost prop- 
agators which have simultaneously the perturbative logarithmic correction and the non-perturbative l/p^ power 
correction. By including the perturbative corrections, the power correction reproduces exactly the leading OPE 
result suggesting the existence of dimension two condensate. 



1. Introduction 

Various condensates are related to non- 
perturbative properties of QCD. The importance 
of the condensates such as (FF) and (qq) are well- 
known. 

Though the dimension two operators such as 
mass terms are not BRST invariant in Yang-Mills 
(YM) theory, it was pointed out recently that 
these operators contain a gauge invariant phys- 
ical part and that a special combination of them 
can be (on-shell) BRST invariant H] [2] [Hi- 
lt was argued that the minimum of along 
the gauge orbit can have a definite physical mean- 
ing, and how to define the physical part non- 
perturbatively . 

There is a BRST-invariant composite operator 
of mass dimension 2 as a linear combination of 

and quadratic ghost, averaged in spacetime: 

-Ai'ix)A^{x) + \iC(x)C(x) 



d X tr 



(1) 



In the Lorentz gauge, it reduces to the same form 
as the known Curci- Ferrari mass term. But we do 
not include it in the Lagrangian. Therefore our 
theory is usual YM. 

It is easily seen that this operator is on- 
shell BRST invariant: SbrstO = 0. The 
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part should be divided into the physical (trans- 
verse) part and other unphysical (longitudinal 
and scalar) part. When we take the vacuum ex- 
pectation value, the CC part precisely cancels 
this unphysical part. Thus the operator has a 
gauge invariant expectation value, though the re- 
maining physical part is nonlocal. Especially in 
the Landau gauge A ^ 0, (O) ^ (iv4^(0)A^(0)). 

Since these operators are not BRST invariant 
as local polynomials, they do not appear in OPE 
of usual gauge invariant quantities. But they may 
appear in OPE of BRST non-invariant quantities 
such as propagators. 

In fact, this operator in the Landau gauge has 
been estimated in various methods. Boucaud et 
al.^ have simulated the lattice propagator and 
used OPE fit to obtain (A^) ~ (lAGeV)^. Ver- 
schelde et al. 0] discussed the effective potential 
of A^. 

Here we take the Schwinger-Dyson (SD) equa- 
tion approach. Recent investigations OIEI of Eu- 
clidean pure SU(A^c) Yang-Mills theory in the 
Lorentz gauge show 



lim p^Dt{p^) 



0, 



lim p^Ggh{p^) = oo. 



(2) 



The transverse gluon propagator vanishes, while 
the FP ghost propagator is enhanced in the in- 
frared limit — > 0. That is, the ghost propaga- 
tor behaves more singularly than the free prop- 
agator in low energy region. These results mean 
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the IR ghost dominance. This is consistent with 
the well-known Gribov prediction and the con- 
finement criterion due to Kugo and Ojima. 

But these results are strongly dependent on ap- 
proximations in the calculation. Simple bare ver- 
tex approximation fails to reproduce 1-loop per- 
turbative result. This situation is very different 
from the case of fermion SD equation. 

In addition to the approximation ambiguity, we 
pay attention to the l/p^ corrections in the UV 
propagators which may come from OPE {{A^) 
condensate) or renormalon or some other non- 
perturbative effects. 

2. SD equations of Yang-Mills theory in 
the Landau gauge 

We consider the SD equation for the gluon 
propagator D and ghost propagator A in pure 
YM theory in the Landau gauge. For the gluon 
propagator D, we adopt the Brown-Penington 



4^ 



to remove a 



projection R^i,{p) := 
tadpole graph. 

Next, we neglect all 2-loop diagrams and adopt 
Higashijima-Miransky approximation to the in- 
ternal propagators. 

Thus the relevant SD equations become 

->-©■■■>•■- = .-y-».-^r.l + ->/'^'^'^>. 
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3. Vertex corrections 



At this stage, in addition to propagators D and 
G, the vertices Fc and F3 are unknown. In the 
Landau gauge, the ghost-gluon vertex Fc does not 
have divergence. So our main attention directs to 
approximate F3. 

In the investigations |S] |H] [7] HI, various 
ansatz are adopted to F3. In [HI, the bare 
or the Slavnov- Taylor (ST) improved vertex is 
adopted. This procedure is similar to the case 
of fermion SD. In [7] |H1j improved vertices are 



used to reproduce the UV 1-loop perturbative re- 
sult correctly where free parameters a and b are 
introduced in [Jj. It does not change UV lead- 
ing exponent. Bloch used a similar ansatz de- 
rived from his own argument based on multiplica- 
tive renormalizability. But in their cases, infrared 
solutions become very delicate in the sense that 
whether ghost dominance property is realized or 
not depends on the choice of parameters or the 
effect of 2-loop diagrams. 

4. UV asymptotic solutions of the SD 
equations 

The simple bare vertex approximation cannot 
reproduce even the 1-loop perturbative result, 
while the IR solutions are very sensitive to the 
choice of vertex correction. In order to get rid of 
this difficulty, we consider the higher logarithmic 
terms in the UV asymptotic solution. We are 
interested in the "non-perturbative" power cor- 
rection terms too. 

To obtain the solutions, we define the gluon 
form factor F and the ghost one G by multiplying 

to the propagators. 

For the UV asymptotic solutions, we adopt the 
new ansatz which has logarithmic (perturbative) 
powers and power (non-perturbative) corrections 
as [5] 



N 



1 ^„ , (lnp2)7+71 ^^^^ 



F(p2) = A(a;lnp2)^^c„(— ^) , ^ 

JV 



G(p2) = B(<^l„p2)« ^ 

^ — ui In 



We use the two parameter ansatz [7] for the 
renormalized triple gluon vertex 



G(g2)i 



^G{(p-g)2)i 



F(q2)i+air((p_g)2)i+b 



(4) 



Simply introducing parameters makes a result 
more ambiguous. In our approach, new parame- 
ters will be automatically determined by consid- 
ering subleading solutions and OPE consistency. 

We put these into SD equations, and get in- 
finite series of algebraic equations with respect 
to coefficients A, B, to, 7, S, Cn, dn, a^^\ fe*^^-*, 71, i5i 
and parameters a, b. 
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5. Result 

We first find tliat without power corrections 

(a^^) = 0, b^^^ = 0) tlie ansatz cannot satisfy 

the coupled SD equation for large , even if we 

include the logarithmic corrections c„, d„ ^ 0. 
A self-consistent solution is obtained when 

7 = -13/22, 5 = -9/44 (5) 

= ^N.XZiAB^ = ^N,X = /3o (A := -|^) 
3 3 167r^ 

for arbitrary value of a and b. 

These reproduce leading exponents and beta 
function in agreement with 1-loop perturbation. 

The remaining coefficients are determined suc- 
cessively except for one degrees of freedom (a + b). 



(a + b) 


421 9152^ 
~ "126 945 ^' 




7 + 71 


9 

= -1.28 -— (a + b), 


(5i = 71 - 1> 


6(1) 


8 




Cl 


35 , 26 , 

= — NcX^ di, 

24 9 






= (4202070705 A^cA^ - 


9583664647VcA2di 




+6683578880AfcAdf 


- 39086260224df ) 




/(304128(645AfcA - 


43648di)), 



C2 = 



: (6) 

Coefficients A and B (overall normalization) will 
be determined if we fix the renormalization con- 
dition, a^^-' and 6^^-' (coefficients of power cor- 
rection) are not determined. Except for di (or 
a + b), all the other coefficients and exponents 
are calculated by simple algebraic equations up 
to any finite orders. The logarithmic expansions 
seem to converge in UV region. 

6. Comparison with OPE result 

The additional logarithmic exponent of the 
power correction calculated from OPE is — (1 — 
^) = —35/44 — —0.795. Corresponding expo- 
nent from the SD solutions in bare vertex approx- 
imation is 7 -I- 7i = —0.935 9 . 

In the improved vertex case only one un- 
known parameter (a -f b) governs this exponent. 
We use the relation ® to determine the param- 
eter. When we set a + b = —4.778 (or di = 



—0.148), the exponent 7 -I- 71 becomes —0.795 
which reproduces the leading OPE. Moreover, 
this value for a + b leads to the ghost dominance 
solution in the IR region. 

7. Conclusion 

We have found the asymptotic SD solutions of 
the gluon and ghost propagators in SU(Af) Yang- 
Mills in Landau gauge for the ansatz using power 
and logarithmic expansion. 

"Non-perturbative" 0(l/p^) power corrections 
are necessary. This corresponds to dimension two 
vacuum condensates [1-3]. 

Contrary to the fermion SD case, the bare ver- 
tex approximation cannot reproduce the UV solu- 
tion consistent with the UV perturbation and the 
UV OPE. We have found the consistent solution 
by considering the UV higher order corrections. 

Using the same vertex ansatz (4) determined 
as above, we have obtained also the IR solu- 
tions which have the same form as the power- 
law solutions [6]: F{p^) ~ A ■ (p2)i-54^ G(j)^) ~ 
B ■ (p^)"''-''^. Therefore the existence of non- 
perturbative power correction does not influence 
the IR solution. IR ghost dominance is realized, 
and the color confinement criterion of Kugo and 
Ojima is satisfied. 
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